Perturbative and Non-perturbative Lattice Calculations for the Study 
of Parton Distributions 

Stefano Capitani a 

a DESY Zeuthen, John von Neumann-Institut fiir Computing (NIC) 
Platanenallee 6, 15738 Zeuthen, Germany 

We discuss how lattice calculations can be a useful tool for the study of structure functions. Particular emphasis 
is given to the perturbative renormalization of the operators. 



1. INTRODUCTION 

Lattice QCD provides non-perturbative tech- 
niques for the computation of the Mellin moments 
of the structure functions of hadrons from first 
principles, without model assumptions. Many re- 
sults have been obtained in the past years from 
the lattice, and include the calculation of the low- 
est moments of various structure functions of the 
quarks: the unpolarized structure functions, the 
spin-dependent structure functions <?i and g2 , and 
the transversity structure function h±. 

Lattice perturbation theory is essential for the 
renormalization of the relevant operators. Per- 
turbative calculations of renormalization factors 
have been carried out for the lowest three mo- 
ments of all structure functions, in the case of 
Wilson and also of overlap fermions (which do 
not break chiral symmetry). For many operators 
these renormalization factors are also known in 
the improved theory. Perturbative renormaliza- 
tion factors have also been calculated for some 
higher- twist matrix elements (4-quark operators). 

2. STRUCTURE FUNCTIONS 

It is not possible to compute a complete struc- 
ture function directly on the lattice (which is set 
up in Euclidean space). The reason is that the 
structure functions describe the physics close to 
the light cone, and this region of Minkowski space 
shrinks to a point when one goes to Euclidean 
space, where Monte Carlo simulations are per- 

*Presented at RADCOR 2002 and Loops and Legs 2002, 
Kloster Banz (Germany), September 8 to 13, 2002. 



formed. However, on a Euclidean lattice it is 
possible to compute the moments of the structure 
functions, using the operator product expansion: 

/ x n F (i) (x,Q 2 ) ~ CM(^) • (h\0^(n)\h). 
Jo M 

The Wilson coefficients contain the short-distance 
physics, and can be perturbatively computed in 
the continuum. The matrix elements contain the 
long-distance physics, and can computed using 
numerical simulations, supplemented by a lattice 
renormalization of the relevant operators. 

Operators of twist two (twist: dimension minus 
spin) dominate the expansion above. Moments 
of the unpolarized structure functions (which 
give the unpolarized distribution q), of the spin- 
dependent structure functions g\ (which gives 
the helicity distribution Aq) and <?2, and of the 
transversity structure function hi (which gives 
the transversity distribution Sq) are measured by 
towers of hadronic matrix elements: 

(x n ) ~ (h\^D fll ---D^\h) 

((Ax) n ) ~ (h\^ 5l{ll D lil ---D lln} ip\h) 

(x n ) g2 ~ (h\ ^IsH^D^D^ ■■•D lin} ^\h) 

((5x) n ) ~ (h\^ 5 a fi{fll D^---D^\h). 

3. OPERATOR RENORMALIZATION 

To obtain physical continuum matrix elements 
from the lattice, the results of Monte Carlo sim- 
ulations need to be renormalized from the lattice 
to a continuum scheme. To obtain these renor- 
malization factors one has to compute 1-loop ma- 
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trix elements on the lattice as well as in the con- 
tinuum. The lattice operators at tree level, for 
p <C 7r /a, have the same matrix elements as the 
original continuum operators, so that at 1 loop 
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where for the continuum one chooses the MS 
scheme, since the Wilson coefficients are known 
in this scheme. In general R l jf ^ Ri° nt , and the 
1-loop renormalization factors on the lattice and 
in the continuum are different. The anomalous di- 
mensions are equal, and thus for ^ = l/a only a 
finite renormalization connects the two schemes. 
The matching between Monte Carlo numbers and 
the renormalized physical results is then 



(or') = (k 
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Since Lorentz symmetry is broken on the lattice 
(which is invariant under the hypercubic group), 
and for Wilson fermions also chiral symmetry, 
the mixing of operators under renormalization is 
more complicated than in the continuum theory. 
Operators multiplicatively renormalizable in the 
continuum may mix with other operators when 
they are put on the lattice. Sometimes these mix- 
ing coefficients are even power divergent, that is 
they behave like l/a™. In general then the match- 
ing factors R]^ — R1° nt are n °t square matrices. 

4. WILSON FERMIONS 

Most calculations are done using the discretiza- 
tion of the (Euclidean) QCD action proposed 
by Wilson. Wilson fermions however break chi- 
ral symmetry, and this causes an additive mass 
renormalization even for a zero bare mass, and a 
heavy pion (m T ~ 500 MeV). Extrapolations to 
the chiral limit are then needed, and this source 
of systematic errors has to be controlled. 

Lattice perturbation theory is quite cumber- 
some: there are more vertices and more diagrams 
than in the continuum, the expressions contain a 
huge number of terms, and the integrals are more 



complicated. The Wilson quark propagator is 



4 E M 7^ sin akp + am f + 2r J2p sin : 
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2r E M sm 



amf 



J2p sin 2 akfi 
which in the continuum limit gives 

5ab -i£^A + "»/ 

E M k*+m 2 f 
and the gluon propagator (in covariant gauge) is 
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The lattice theory has exact gauge invariance 
at any finite a. This causes the presence in the 
action of the group elements Up — exp(iag Ap) 
instead of the algebra elements Ap . One has then 
to expand the UpS in terms of the Ap's, and this 
generates an infinite number of vertices, of which 
only a finite number is needed at any given order 
in go. The quark-quark-gluon vertex is 
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(t a ) bc (i 



n> co s — + r sin 
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(where s = p\ + P2, the momenta of the quarks), 
which for a — » gives — go(t a ) bc i / y fl , while the 
quark-quark-gluon-gluon vertex is 

-\agl V{t», t b V d (~h, sin^+r cos ^f) .(2) 

This is an irrelevant vertex (is zero in the contin- 
uum limit), but still gives non- vanishing contri- 
butions to Feynman diagrams in divergent loops. 

Also in the pure gauge part there is an infinite 
number of vertices. The 3-gluon vertex is 

• a(k-p) p aq^ 



i5o/ ab M<Wsm 

a 2 



cos ■ 



. a(p - q)p ak v 
+5 vp sin — £ cos — 



+5 0U sin 



. a(q - k) v ap p 



cos • 



which in the continuum limit becomes 
i.9o/ abc {<W fc - p)p + Mp - l)n + <W(? - k )»}- 
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The 4-gluon vertex is too complicated to be re- 
ported here. Furthermore, the gauge measure at 
order <?q gives a 1/a 2 mass counterterm. Finally, 
using a Faddeev-Popov procedure one can obtain 
the Feynman rules for the ghost propagator and 
the ghost interactions. The effective ghost-gauge 
field interaction is not linear in the gauge poten- 
tial A^, and thus also in this sector new vertices 
appear that have no continuum analog, like the 
ghost-ghost-gluon-gluon vertex. 

5. OVERLAP FERMIONS 

A Dirac operator D — J^D^ which satisfies the 
Ginsparg- Wilson relation J5D+DJ5 = a/pDj^D 
defines fermions with maintain an exact chiral 
symmetry also for non-zero lattice spacing and 
also keep all other fundamental properties. 

One of the possible solutions of the Ginsparg- 
Wilson relation is given by overlap fermions. In 
the massless case the overlap-Dirac operator is 
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X = D 
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where Dw is the Wilson-Dirac operator. For < 
p < 2r one has the right spectrum of massless 
fermions. Additive mass renormalization is then 
forbidden, and one avoids a source of systematic 
errors always present using Wilson fermions. 

The massless quark propagator in the overlap 
is more complicated than in Wilson, and is 



\2p (uj(k) + b(k)) 
where b(k) = 1/a (2r sin 2 afc^/2 — p), 
uj{k) = ^xl(k)X Q (k), and 

X (k) = - (i ^2 sin ak v + 2r^sin 2 ^ - 




fj, fj. 

The overlap vertices can be expressed in terms 
of the vertices of the QED Wilson action Wi^ and 
W 2fi , which are the vertices (Q) and (||) without 
the color matrices, and of the quantity Xq. For 
example, the quark-quark-gluon vertex is 



{t a ) bc ■ p 



1 



x (p 2 )wi(p 1 ,p 2 )x (p 1 : 



The quark-quark-gluon-gluon vertex is very long 
and will not be given here (see for example [ EJ). 

6. IMPROVEMENT 

It is very expensive to decrease the errors due 
to the granularity of the lattice just by reducing 
the lattice spacing a, because the cost would grow 
like a~ 5 in the quenched theory, and in full QCD 
it would grow even faster. Halving the discretiza- 
tion errors in this way would then be at least 30 
times more expensive. A more effective way is to 
improve actions and operators. 

0{a) improvement removes the contributions 
of order a to the systematic error arising from 
the finiteness of the lattice spacing by adding a 
counterterm to the fermion action, so that 
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For Wilson fermions, this is achieved in on-shell 
matrix elements adding a counterterm whose co- 
efficient c sw has to be exactly tuned. Only for its 
appropriate value, for a given go, the 0(a) effects 
are canceled and one gets a faster convergence to 
the continuum limit. Due to this action countert- 
erm one has to add to the Wilson vertex (|l|) the 
improved quark-quark-gluon interaction vertex □ 

-c sw -go^{t) cos ^a^x smaqx, 

A 

where q is the difference between the incoming 
and outgoing momenta of the quarks. In addi- 
tion to this, one also has to improve the various 
operators, that is one must add a basis of higher 
dimensional (irrelevant) operators with the same 
quantum numbers, with dim(Oi) = dim(O) + 1: 



O — ► O imp = (1 + b o am)0 



The 0(a) improvement for overlap fermions is 
much simpler. The action is already improved, 



2 The fermion propagator and the vertices with an even 
number of gluons are instead not modified by the improve- 
ment. Neither is the gluon propagator: the first correc- 
tions to the pure gauge action are already of order a 2 . 
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and thus there is no need of new interactions. The 
improved operators are given, to all orders, by 

0^ = ^1-^)0(1-^)^ 



(while for Wilson fermions the improvement con- 
struction is different for each operator) . One thus 
gets full 0(a) improvement without tuning any 
coefficients, and to all orders of perturbation the- 
ory. For Wilson fermions, instead, one has to 
determine the coefficients of the operator coum- 
terterms order by order in perturbation theory. 13 

7. FORM CODES 

Due to the complexity of the calculations, to 
the possible great number of diagrams, and above 
all to the huge amount of terms in each diagram, 
computer codes have to be used. To evaluate the 
Feynman diagrams and obtain the algebraic ex- 
pressions for the renormalization factors, we have 
thus developed sets of computer codes written 
in the symbolic manipulation language FORM. 
These codes take as input the Feynman rules for 
the particular combination of operators, propa- 
gators (Wilson or overlap) and vertices (Wilson, 
improved, or overlap) appearing in each diagram, 
expand them in the lattice spacing a at the ap- 
propriate order, evaluate the gamma algebra on 
the lattice, and then work out everything until 
the final expressions are obtained. The main dif- 
ficulty in developing such codes is that the (Eu- 
clidean) Lorentz group 0(4) breaks down to the 
hypercubic group W4. A representation which is 
irreducible under the (Euclidean) Lorentz group 
is thus in general reducible under the symmetry 
group of the lattice. This gives rise to a whole 
new class of problems, of which the most serious 
concerns the Einstein summation convention. 

The FORM language has been developed hav- 
ing in mind the usual continuum calculations. 
There are therefore many useful built-in features 
which are sometimes somewhat of an hindrance 
when one tries to perform lattice perturbative cal- 
culations. A blind use of them would give, for 



3 This appears to be a difficult task, and even for the first 
moment of unpolarized distributions all lowest-order im- 
provement coefficients have not yet been determined. 



i sin fc> 



example, 

^ 1\P\ sin k x 
x 

One thus needs to develop special routines to cope 
with the gamma algebra on the lattice. 

Integrals are done over the first Brillouin zone 



d 4 k 



f(ak,ap) 



a 4 f(k',ap) 



with the rescaling k' = afc, and in general an ex- 
pansion in ap is needed. FORM codes become 
then necessary also because of the huge number 
of terms arising from the Taylor expansions of the 
Fourier transforms of operators, propagators, ver- 
tices and covariant derivatives. The n-th moment 
of a parton distribution behaves like 

®YD^.--D^) ~ JL 

because D ~ 1/a, and thus one has to perform an 
expansion in a to order n for every quantity. For 
example, the Wilson quark propagator to 0(a) is 



5 ab D- 1 (k)i (-i^7 M sin^ + 2rJ]sh 
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-D-\k) 
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r E ^ sin kfj 



i 7 P sin k p + 2r 
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where the denominator is 
D{k) = J^sin 2 fc M + (2r 

As a consequence, a huge number of terms ap- 
pear (at least in the initial stages of the manip- 
ulations), and also the gamma algebra becomes 
quite cumbersome to do by hand. All this also 
implies a limitation on the number of moments 
of structure functions that one can compute. □ 



This comes on top of the limitation coming from operator 
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8. RESULTS 

We have calculated the renormalization in the 
MS scheme of several operators measuring the 
lowest moments of structure functions [ Q . In 
some cases for a given moment we have computed 
two operators, labeled (a) and (b), which belong 
to two different representations of the discrete Eu- 
clidean Lorentz group (see Table 0). eI Operator 
mixing for overlap fermions is simpler than for 
Wilson fermions: chiral symmetry is not broken 
and thus it prohibits any mixing with operators of 
different chirality. Mixing coefficients which are 
power-divergent like a~ n in the continuum limit 
can be eliminated from the start from overlap cal- 
culations if the corresponding operators belong 
to multiplets with the wrong chirality. In the 
Wilson case, indeed, the operators measuring the 
moments of the gi structure function present ad- 
ditional mixings with wrong-chirality operators 
with power-divergent coefficients. Furthermore, 
the Z'b of the n-th moments of q and Aq are not 
constrained to be equal anymore. 

For reasons of computing power, structure 
functions have been studied mostly in the 
quenched approximation, and only recently com- 
putations in full QCD have become feasible. On 
the lattice the Grassmann variables have to be an- 
alytically integrated, and the partition function 




becomes in the simulations 

Z = JvUdet(p[U] + mf)e~ s " lu] , 

which is equivalent to use an effective action 
S eS [U] = S g [U] - lndet(#[C/] +m f ). 

Quenching amounts to doing simulations with 
det(p[U] + 77i/) = 1. This means that there are 
no sea quarks: all internal quark loops are ne- 
glected. Although it looks quite drastic, it is not 
a bad approximation for many physical quanti- 
ties, and apparently (as it has been understood 
recently) also for structure functions. 

mixing: when an operator has more than four indices, two 
of them have to be equal, and this opens the door to lots 
of mixings, often with power-divergent coefficients. 
5 In the operators (a) all Lorentz indices are distinct. 
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Table 1 

Renormalization constants of multiplicatively 
renormalized structure function operators, for 
a = 0.1 /m, in the MS scheme. 



Structure functions in full Wilson QCD have 
been the major progress in the past couple of 
years [ |[ The first studies with dynamical 
Wilson quarks have shown no statistically signif- 
icant differences between quenched and full QCD 
results. Previously it had been conjectured that 
quenching was the cause of the observed discrep- 
ancies of lattice results with experiment. Evi- 
dently this is not the case, at least for the values 
of Wilson quark masses attainable at present. 

The largest discrepancies appear for unpolar- 
ized quark distributions. The spin-dependent g\ 
and g2 structure functions have also been exten- 
sively studied, as well as the axial charge qa = 
Ait — Ad (which deviates less from experiment) 
and the hi transversity structure function [ 
whose lowest moment is the tensor charge 8u—8d. 

One limitation of present simulations is given 
by the difficulty to compute numerically diagrams 
which are disconnected (except through gluon 
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lines). The disconnected diagrams are however 
flavor independent, and do not contribute to the 
difference between u and d structure functions 
(for SU(2) degeneracy). This means that quanti- 
ties like gA = Ait— Ad and {x n ) u -d do not receive 
contributions from disconnected diagrams. Their 
lattice results seem also to be closer to the exper- 
imental numbers. 

The chiral extrapolations of the lattice data lin- 
ear in m q ~ rrv^ could also be responsible for the 
discrepancies with experiment. It has been re- 
cently suggested that extrapolations using chiral 
perturbation theory could solve these discrepan- 
cies (see [ D and references therein). In fact, the 
pion cloud of the nucleon is not adequately de- 
scribed by current lattices. Extrapolation formu- 
lae which use chiral perturbation theory look like 
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where A is a phenomenological cutoff related to 
the size of the source generating the pion cloud. 

These extrapolation formulae seem able to re- 
solve the discrepancy with experiment. However, 
results can be reproduced only with A's with a 
range so wide to have no predictive power, and 
the currently available lattice data do not even 
discriminate between linear and chiral perturba- 
tion theory fits. The problems is that the avail- 
able pions are not sufficiently light. A smaller 
pion mass (m w < 250 M eV) is needed before the 
parameters of the chiral expansions can be well 
determined on the lattice. Such calculations re- 
quire about 8 Teraflops for one year. The next 
generation of computers, coming in a couple of 
years, will be able to perform these calculations. 

There is no doubt that the pion cloud of the 
proton is very important. For it to be adequately 
included in the lattice box and properly measured 
in Monte Carlo simulations, the pion correlation 
length should be much smaller than the lattice 
size. This points to the use of larger lattices. 

Higher-twist corrections have also been stud- 
ied [ 0. The lattice results for the twist four 
(l/Q 2 ) corrections are (in a few particular cases) 
quite small. The renormalization factors of four- 
quark operators like ^ a tp^^^t""^ ■ / ipji / j5t a ip 



have been calculated in perturbation theory, and 
the twist-4 results for the first moment of the 
unpolarized pion and proton structure functions 
turn out to be much smaller than the correspond- 
ing twist-2 matrix elements computed on the lat- 
tice, and also smaller than the phenomenologi- 
cal numbers. These are not however complete 
and systematic studies, like the ones concerning 
leading-twist operators. Only particular flavor 
and isospin combinations could be considered, to 
avoid mixing with lower-dimensional operators. 

9. CONCLUSIONS 

The lattice provides invaluable techniques for 
investigating moments of structure functions non- 
perturbatively from first-principles. The match- 
ing to the MS scheme is done by computing renor- 
malization factors in perturbation theory (Wilson 
and overlap). Simulations with Wilson fermions 
are now also performed in full QCD, but many 
discrepancies between lattice and experiment re- 
main. More control on the systematic errors is 
required. 
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